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Abstract 

In a R-parity violating (RPV) model of neutrino mass with three bilinear cou- 
plings Hi and three trilinear couplings A^ 33 , where i is the lepton index, we find 
six generic scenarios each with a distinctive pattern of the trilinear couplings 
consistent with the oscillation data. These patterns may be reflected in direct 
RPV decays of the lighter top squark or in the RPV decays of the lightest su- 
perparticle, assumed to be the lightest neutralino. Typical signal sizes at the 
Tevatron RUN II and the LHC have been estimated and the results turn out to 
be encouraging. 
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1. Introduction 



Neutrino oscillations have been observed in different experiments [1], from which it is 
confirmed that the neutrinos have tiny masses, several orders of magnitude smaller than any 
other fermion mass in the Standard Model ( SM ). If right handed neutrinos are introduced 
in the SM one can formally get Dirac masses of the neutrinos. But the corresponding Yukawa 
couplings must be unnaturally small. 

The SM being unable to naturally explain the origin of very small neutrino masses , Im- 
parity violating (RPV) supersymmetry (SUSY) [2, 3] could be a viable alternative. In the 
lepton number violating version of these models, the neutrino masses (Majorana type) and 
mixing angles [4, 5] can naturally arise without requiring the right handed neutrinos. 

In the minimal supersymmetric extension of the SM (MSSM) there is also a R-parity 
conserving (RPC) sector with Yukawa like couplings and soft breaking masses/ mass param- 
eters, which govern the masses of the superpartners of the SM particles (the sparticles). As 
demanded by the naturalness argument, these masses should be of the order of one TeV. 
The theoretical predictions for the neutrino masses and mixing angles depend on these RPC 
parameters, in addition to the parameters of the RPV sector. Thus the exciting program of 
sparticle searches and the reconstruction of their masses at the on going ( Tevatron RUN II) 
and the upcoming (the Large Hadron Collider (LHC) or the International Linear Collider 
(ILC)) accelerator experiments have the potential of testing the RPV models of v- mass. 

In R-parity conserving( RPC) SUSY the lightest supersymmetric particle (LSP) is pre- 
dicted to be stable and a carrier of missing transverse energy [fir). In contrast the model 
under consideration allows the LSP to decay via RPV couplings into lepton number violating 
modes. The multiplicity of particles in any event is, therefore, much larger compared to the 
corresponding event in RPC SUSY containing the stable LSP leading to distinct collider 
signatures. Moreover, since the LSP is not necessarily a carrier of fir the reconstruction of 
sparticle masses appears to be less problematic than that in the SM. 

Another characteristic signature the model under consideration is the direct decays of 
sparticles other than the LSP into lepton number violating channels. Apparently the strin- 
gent constraints on the RPV couplings obtained from the neutrino data (see below) suggest 
that the branching ratios(BRs) of these decays will be highly suppressed compared to the 
competing RPC decays. One notable exception is the direct RPV decay of the lighter top 
squark [6, 7, 8, 9], if it happens to be the next to lightest supersymmetric particle(NLSP). 
This assumption is theoretically well-motivated due to large mixing effects in the top squark 
mass matrix. In this case the RPC decays of the top squark are also suppressed (to be 
elaborated later) and can naturally compete with the RPV decays even if the underlying 
couplings are as small as that required by the neutrino data. Thus the competition among 
different decay modes of the lighter top squark, which may be observed during Tevatron 
RUN II, is a hallmark of RPV models of neutrino mass [9] . 

Of course the see saw mechanism [10] which can be naturally implemented in any grand 
unified theory(GUT) not necessarily supersymmetric, provides an elegant alternative expla- 
nation of small neutrino masses. Unfortunately the simplest version of this theory, a GUT 



1 



with a grand desert, predicts a low energy spectrum practically identical with that of the 
SM. Thus there is no testable prediction outside the neutrino sector. 

In view of the large number of RPV parameters, constraints on them [3, 11] are usually 
obtained from the experimental data by employing a simplifying assumption. It is assumed 
that only a minimal set of the parameters contributing to the observables under study 
are numerically significant. Following this approach one usually analyzes some benchmark 
scenarios, each consisting of a minimal set of RPV parameters at the weak scale [12, 13], 
capable of reproducing a phenomenologically viable neutrino mass matrix. Constraints on 
the parameters belonging to each scenario are then obtained by using the neutrino oscillation 
data [12, 13]. 

Among the examples in Ref. [12], we have focused on a specific model with three trilinear 
couplings A^ 33 ( where i = 1,2,3 is the lepton generation index ) and three bilinear RPV 
parameters /ij.The stringent upper bounds on the trilinear (bilinear) couplings are ~ 1CT 4 
( ~ 10~ A GeV ). As a result the contributions of these couplings to most low energy processes 
except LSP decays are negligible. As already mentioned, a notable exception could be the 
direct RPV decay of the lighter top squark [8, 9] into a 6-quark and a charged lepton. 
Moreover, by reconstructing the lepton - jet invariant mass the lepton number violating 
nature of the decay can be directly established [9]. A model independent estimate of the 
minimum observable branching ratio (MOBR) of the channel t\ — > e + d as a function of the 
lighter top squark mass (m^ ) for Tevatron RUN II was also presented [9] . This estimate was 
then translated into the magnitudes of the A's for representative values of the parameters in 
the RPC sector and they were found to be close to the bounds obtained from the oscillation 
data. 

In this paper we go a step beyond simple estimates and obtain testable quantitative 
predictions. We find various combinations of three A'-type couplings allowed by the neutrino 
data. We do this by randomly generating 10 9 sets of the above six RPV couplings for 
representative values of the parameters in the RPC sector and identify the sets consistent 
with the oscillation data [14]. Our numerical results are checked by analytical calculations 
in a simple approximation (see sections 2 and 3 for the details). An interesting common 
feature of all the allowed sets is that one fii is much smaller than the other two and for 
each such scenario there are two characteristic patterns of the A^ 33 couplings. We thus have 
six generic scenarios. In each scenario the relative BRs of the three decay modes t\ — > l^b 
reflect the underlying model. Thus if some of the RPV decay channels of ti are observed with 
BRs as predicted, a particular generic scenario can be vindicated. Of course more definite 
conclusion can be drawn if complementary information about the RPC sector, the masses of 
the sparticles in particular, can be obtained kinematically (by measuring invariant masses 
of the final state particles, edges of the energy distributions of some of the decay products 
etc). Some simple examples will be given later. 

In our numerical work, we have chosen the parameters of the RPC sector corresponding 
to several scenarios classified according to the properties of the electroweak gaugino sector 
which determine the tree level neutrino masses ( see section 4 for the details). In each case the 
lightest neutralino is assumed to be the LSP. We then study the predicted competition among 
various RPC [15, 16] and RPV decay modes of the top squark NLSP for RPV couplings 
allowed by the neutrino data. We find that many of the solutions predict BRs larger than 
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the estimated MOBR at Tevatron in [9] . We have also calculated the number of some typical 
signal events, where the expected background is low from top squark production and decay 
at Tevatron RUN II and LHC. 

If on the other hand no signal is seen during Tevatron RUN II the allowed parameter 
space(APS) will be significantly squeezed. From the limits on the trilinear couplings obtained 
from Tevatron RUN I and RUN II [8] and from the projected sensitivity of RUN II data to 
these couplings [9], it seems that A- 33 ~ 10~ 4 can be ruled out for lighter top squark masses 
within the kinematic reach of the Tevatron. 

If the top squark is not the NLSP and the RPV couplings are as small as that required 
by the current neutrino data [14], the decay of the LSP would be the only signature of the 
RPV model of neutrino mass. However, the amount of information that may be extracted 
will depend on parameters of the RPC sector, in particular on the LSP mass. If the LSP is 
lighter than the top quark, then it decays via the modes 

x? -> nbb (l) 

where I — e , /j, or r. Since the neutrinos are not detectable the branching ratios of the 
individual channels can not be measured. The decay length of the LSP depends on many 
model parameters belonging to the RPC and RPV sectors and pin pointing the underlying 
model of neutrino mass from one observable may not be an easy task (see, however, [17]). In 
addition the lepton number violating nature of the underlying interaction can not be directly 
established since the neutrinos escape the detector. 

The signatures will be unambiguous if the Xi happens to be heavier than the top quark. 
In this case apart from the decay channel in Eq. (1) the LSP will also decay into a charged 
lepton ( e, /i or t ), the top quark and the anti-bottom quark, a clearly lepton number 
violating signature. This observation motivates us to also scan the parameter space where 
the LSP is heavier than the top quark. Of course the decay modes involving the top quark will 
be phase space suppressed compared to the ones in Eq. (1). Nevertheless our computations 
show that the branching ratios of the three modes involving the t-quark are numerically 
significant over the entire parameter space allowed by the oscillation constraints. Since the 
LSP is rather heavy in this scenario it is unlikely to be produced at the Tevatron. However, 
observation of all four modes and measurement of their branching ratios at the LHC or ILC 
will provide crucial tests of the underlying model of neutrino mass. We have calculated the 
size of some typical signals from the production of electroweak gaugino pairs at the LHC. 

The plan of the paper is as follows. In section 2 we establish the notations and briefly 
review the neutrino mass matrix in the RPV model under study. In section 3 we identify the 
six generic scenarios of the RPV sector compatible with the neutrino data and using several 
representative values of the parameters of the RPC sector obtain sets of RPV parameters 
allowed by the oscillation data. In section 4 the top squark decays are studied in different 
scenarios. The LSP decays are analyzed in section 5. Our conclusions and future outlooks 
are summarized in section 6. 
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2. Neutrino mass matrix at the tree and loop level 



R-parity is a multiplicative quantum number defined as follows [3], 

R p = (_l)3B+L+25 ; (2) 

where B is the baryon number, L is the lepton-number and S denotes the spin. For particles 
R p — +1 and for sparticles R p = -1. 

In general the superpotential of the MSSM may contain RPV terms which violate both 
B and L conservation. This leads to catastrophic proton decay with a mean life time not 
allowed experimentally. All B and L violating terms can be removed from the superpotential 
by imposing R-parity as a symmetry. The resulting model is known as the R-parity con- 
serving MSSM. In order to prevent proton decay it is, however, sufficient to remove either 
B-violating or L-violating terms by imposing appropriate discrete symmetries. Models with 
B-violating terms only can not generate neutrino masses. As discussed in the introduction 
we have focused on a specific model of neutrino mass with A' type couplings only. The 
general R-parity violating superpotential of our interest takes the form: 



W = Wmssm + W^, (3) 
W K = K 3k e ab LtQ b 3 D c k + fi t e ab L^H b 2 . (4) 

Here, Wmssm[2] is the usual superpotential of the MSSM containing the terms which give 
mass to the SM fermions. The i, j, k — 1,2,3 are generation indices and a, b = 1,2 are SU(2) 
indices and V denotes charge conjugation. The A's are dimensionless trilinear RPV Yukawa 
like couplings, /ij's are bilinear RPV terms [3, 18] with dimensions of mass, which determines 
the amount of mixing between the lepton and Higgs superfields. In Eq.(4) L, Q and Hi 
denote, respectively , SU(2)l doublet lepton, quark and up type higgs superfields and D is 
the SU{2) L singlet down type quark superfields. One can also construct models of z/-mass 
with RPV A-type couplings [3, 4, 5, 19] only. The phenomenology of such models will be 
quite distinct from that discussed in this paper. 

The tree level and loop level neutrino mass matrices are given below [3, 5, 12]: 



MX C = C W , (5) 

and C is given by : 

2 {M 1 + tWL 2 6 w M 2 ) 2 

C = g > IdelM Vd (6) 

where gi is the SU{2) gauge coupling, Mi,M 2 are the U(l) and SU(2) gaugino masses 
respectively and det M is the determinant of the R-parity conserving neutralino mass matrix. 
Here we are working in the basis where the sneutrino vacuum expectation values (vevs) are 



| + vl, where (v u ) is the vev of down (up) type Higgs field. The one loop 
mass matrix is given by : 
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M loo P 



K 2 Xf 33 

-^2A 13 3A 2 33 
-^2A 13 3A 3 33 



K 2 X' 133 X'233 

K 2 X 233 
K 2 X 233 X 333 



^2X' 133 X' 333 

-K2A233A333 

K 2 X 
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(7) 



2^333 



where K 2 is given by : 



, X b f{Xg) 

' 16vr 2 M { q f 



R) . 



(8) 



with 



/(*) 



In x 



x 



M, 



(3) 



and 



(9) 



X b = A 6 -//tan/3. (10) 

In the above X b is the off-diagonal mixing term in the b-squark mass matrix , A b is the 
soft trilinear term for the bottom squarks, tan/? is the ratio of v u to Vj, /i is the higgsino 
mass parameter , m b is mass of the bottom quark . M^ q2 are the masses of the two b - 
squark mass eigenstates. 

In this work we shall assume that the masses of the right handed and left handed squarks 
are equal. The same assumption applies to the slepton sector. 
Thus, we take mass matrix up to the one-loop to be 

M u = MT C + M l ° op . (11) 

It has been noted in the literature that there may be other loop contributions to the 
neutrino mass matrix [20, 5]. For example, the soft breaking RPV bilinear terms may 
contribute to some of the loops [5]. As has already been mentioned the large number of 
free parameters compels one to work in benchmark scenarios with a limited number of RPV 
parameters. 

In Table [1] we present the neutrino data that has been used for the numerical work 
in this paper [14]. We shall consider the data at the 2a level. The notations used are 
as follows [21]: the neutrino mass squared differences are A 2 sol = Am 21 = \m\~m\\ and 
Am 2 tm = Awijj = \m 3 — m\\ respectively, where m 1; m 2 and m 3 are the three eigenvalues 
of the neutrino mass matrix in Eq. 1 [12]. The mixing angles are extracted from the 
eigenvectors corresponding to appropriate mass eigenvalues. 



3. The six generic scenarios allowed by neutrino data 

The neutrino mass matrix in section 2 can be recasted in the following form: 
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parameter 


best fit 


2a 


3a 


Aa 


AmL[10- 5 eF 2 ] 


8.1 


7.5-8.7 


7.2-9.1 


7.0-9.4 


sin 2 9 12 


2.2 


1.7-2.9 


1.4-3.3 


1.1-3.7 


0.30 


0.25-0.34 


0.23-0.38 


0.21-0.41 


sin 2 9 2 3 


0.50 


0.38-0.64 


0.34-0.68 


0.30-0.72 


sin 2 #13 


0.000 


< 0.028 


< 0.047 


< 0.068 



Table 1: Best-fit values, 2a, 3a, and Aa intervals for the three flavour neutrino oscillation 
parameters from global data analysis [14] including solar, atmospheric, reactor (KamLAND 
and CHOOZ) and accelerator ( K2K ) experiments [1]. 



M, 




(12) 



where Tj and Di are given below : 



D i = Cfi 2 i + K 2 Xf 33 (13) 

D 2 = Ctf + K 2 X' 2 \ 3 (14) 

D 3 = C^ + K 2 Xg 3 (15) 

7\ = C/ix/12 + K 2 X[ 33 X' 233 (16) 

T 2 = Cfi m + K 2 X' 133 X' 333 (17) 

T 3 = Cfi 2 fi 3 + K 2 X' 233 X 333 (18) 



The oscillation parameters can be easily calculated analytically if any one of the TjS vanish 
and Di ^ for all i. Thus we tentatively propose the following generic scenarios : 

(a) nx = and either A^ 33 = (7\ = 0) or X' 333 = (T 2 = 0) 

(b) H2 = and either A' 133 = (7\ = 0) or A 333 = (T 3 = 0) 

(c) fi 3 = and either A' 133 = (T 2 = 0)or A' 233 = (T 3 = 0) 

The above patterns help us to classify different regions of the RPV parameter space 
consistent with the oscillation data in Table [1]. 

For illustration we consider a particular hierarchy ^ <^ fi 2 or fi 3 (generic scenario 
(a) ) and the following representative choice of RPC parameters : Mi=110, M2=200, // 
=400, tan/9 =5, M^=400, A^=1000 where all masses and mass parameters are in GeV. We 
randomly vary all six RPV parameters within the ranges shown in Table [2] columns 2 and 
3 and generate 10 9 sets of parameters. Only jii is constrained to be rather small. 

We then pick up one solution corresponding to ji\ = 0.08 x 10 4 eV, /i 2 = 140 x lO 4 ^, /i 3 = 
111 x 10 4 eV, X' 133 = 3.94 x 10~ 5 , X' 233 = .0 1 5 x 10 ~ 5 and A' 333 = 7.03 x 10" 5 . The numerically 
calculated oscillation parameters are presented in the second column of Table [3] . In column 
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Given 


ranges 


Allowed ranges 


RPV parameters 


Max 


Min 


Max 


Min 


Hi[10 4 eV] 
^ 2 [l0 4 eV] 
/i 3 [10 4 eV] 

a; 33 [io- 5 ] 

A 333 [10- 5 ] 


50.0 
300.0 
300.0 
12.0 
12.0 
12.0 


0.001 
1.0 
1.0 
0.001 
0.001 
0.001 


19.0 
134.0 
178.0 
4.32 
9.52 
9.67 


0.10 
90.0 
89.0 
3.2 
0.01 
0.002 



Table 2: Allowed ranges of RPV parameters in scenario (a) 

3 of the same table we present the analytically calculated results with the approximation 
Hi = 0, A 233 = 0, while other RPV parameters are as above. The agreement between the 
two sets provides a test of the reliability of the numerical procedure. Changing the range of 
variation of /ii does not lead to any qualitatively new solution with /ii comparable to /i 2 or 
^3. Thus classifying a generic scenario by a a small magnitude of ji\ is indeed valid. 



Neutrino oscillation 


Numerical results 


Analytical results 


parameters 






Am 2 sol [W- b eV 2 } 


8.65 


8.80 


AmLJlO-W 2 ] 
sin 2 9 12 


1.74 


1.69 


0.329 


0.327 


sin 2 9 23 


0.437 


0.435 


sin 2 0i3 


0.016 


0.016 



Table 3: Comparison of Numerical and Analytical results 

Out of 10 9 sets of randomly generated RPV parameters only a few satisfy the data in 
Table [1]. This illustrates the highly restrictive power of the currently available oscillation 
data inspite of the relatively large errors. The allowed ranges of the six RPV parameters 
are presented in Table [2] columns 4 and 5. Examining the entire APS corresponding to 
the hierarchy jii <C /i 2 , /i 3 , where /z 2 , /i 3 can be comparable, one can identify two subclasses 
under the generic scenario (a). The trilinear RPV couplings in the two subclasses are found 
to satisfy the following patterns: 

( a i) A333 > A' 133 > A 233 , 

{ a 2> A 23 3 > A' 133 > A333. 

In fact most of the APS a\ (a 2 ) corresponds to A 233 <C A' 133 ( A 333 A' 133 ). However, 
there are exceptions albeit for relatively small regions of the APS, where the two smaller 
couplings could be of comparable magnitude. 



7 



For ji 2 <C Hi or /i 3 (scenario b) and /i 3 <C /ii or /i 2 (scenario c) each region of the APS 
is also characterized by a specific pattern of the trilinear couplings due to the constraints 
imposed by the neutrino data. If some direct RPV decay modes are observed, these patterns 
would be reflected in the measured BRs revealing the model underlying neutrino oscillations. 

For ji 2 <C Hi or /i 3 , the trilinear couplings follow the patterns 

^133 ~ -^233 >> A333, 

(M A 233 > A333 >> A' 133 . 

For /i 3 <C Hi or ji2, on the other hand, we have 

( c i) A333 > A' 133 >> A 233 , 
(c 2 ) A'333 > A 233 >> A' 133 . 

It can be readily checked analytically that even in the most general case when none of 
the RPV parameters vanish one eigenvalue is still zero. Thus analytical solutions are still 
possible. The formulae for the masses and the mixing angles are somewhat cumbersome. 
For scanning the parameter space we therefore prefer the numerical method. 

It should be noted that there is no straightforward way of determining the Hi parameters 
directly from collider signals. Thus it is gratifying to note that the scenarios {a{) - (c 2 ) can 
be identified by the decay branching ratios alone provided RPV decays into charged leptons 
are observed. 

On the face of it the scenarios (ai) and (ci) look similar. But scanning the entire APS 
in both cases we have found that 

2.6 % < BR(e) < 4.2 % and 19.7 % < BR(r) < 10.2 % in (ai), 

while 

9.5 % < BR(e) < 12.8 % and 18.4 % < BR(t) < 22.0 % in (ci), 

where BR(e) (BR(t)) refer to the BR of any direct RPV decay mode into a final state with 
e (r). Thus each scenario will have its characteristic decay pattern. 

4. The lighter top squark decay 

We now study the collider signals that may be triggered by the sets of A' couplings consistent 
with the neutrino data[14]. The signals can be classified into a few patterns corresponding 
to the six generic scenarios discussed in section 3 and the hierarchy of trilinear couplings 
associated with them. All allowed sets would lead to the RPV decays of the lighter top 
squark (tx) in Eq. (19) with appreciable branching ratios if it happens to be the NLSP, 
which is the case over a large region of the RPC parameter space. 
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(a) h^l+b (19) 

where i = 1-3. From section 3 it is clear that the relative BRs of the three leptonic modes 
will be different in generic scenarios al) - c2). Hence identification of the RPV parameters 
underlying the model of neutrino mass in future experiments is a distinct possibility. 
The RPC decay modes of the lighter top- squark (ti) are listed below : 



(b) h^bxt (20) 

(c) h -> Wu, btv, bWxi (21) 

(d) h^cxl (22) 

(e) h^bxiff (23) 



/ and /' being a quark- ant iquark or l-ui pair. If the lighter top squark is the NLSP, only the 
decay modes d) and e) and the last mode of c) are allowed. The mode in c) will be phase 
space suppressed for t\ masses within the kinematic reach of the Tevatron, which are the 
main subject of this study. Moreover, it will be highly suppressed if the LSP happens to be 
bino like, which is quite natural in popular models like the ones with a unified gaugino mass. 
We shall not consider this mode in this paper. The channels in d) and e) have naturally 
suppressed widths and can very well compete with each other [16] or with the RPV mode, 
especially if A^ 33 is ~ 10~ 4 — 1CT 5 as required by neutrino data [8, 9]. 

Our choices of the RPC parameters are guided by the gaugino sector. It is clear from 
Eq.(5) section 2 that the parameter C sets the scale of the tree level neutrino mass matrix. 
This C depends solely on the parameters of the gaugino mass matrices. Accordingly we have 
chosen the following scenarios. 

1. Models in which the lighter chargino (xf) and the two lighter neutralinos (xi an d Xz) 
are higgsino like (Mi,M 2 /i) and all have approximately the same mass ( ~ //). Thus 
it is difficult to accommodate the top squark NLSP without fine adjustments of the 
parameters. Thus the LSP decay seems to be the only viable collider signature which 
will be discussed in the next section. 

2. Models in which xf, Xi an d X2 are gaugino like (Mi < M 2 <C //) and the top squark 
is the NLSP. 

3. Models in which xf an d X2 are mixed (Mi< M 2 ~ //) and the top squark is the NLSP. 

In all models the parameters are so chosen that the lightest neutralino (xi) happens 
to be the LSP. Further the squarks belonging to the right and left sectors of all flavours 
are assumed to be mass degenerate. In fact although the common squark mass, \i and tan 
j3 occur in both the neutrino and the top squark sector, one can choose the soft trilinear 
parameter A t , which does not affect the neutrino sector, to satisfy the top squark NLSP 
criterion in most cases. However, attention must be paid so that large values of A t do not 
lead to a charge color breaking ( CCB ) vacuum [22]. The BRs of top the squark decay 
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modes and competition among the RPV and RPC decays would be highly indicative of the 
underlying model. 

For model 2, the gaugino like model, we use the following representative values of the 
parameters of the RPC sector: Mi = 110, M 2 = 200, fj, = 400, tan/3 = 5, Mg = 400, 
At, = 1000, A t = 970 and Ma = 300 where all masses and mass parameters are in GeV and 
Ma is the CP odd higgs mass. The first six parameters, which are the same as the ones used 
in section 3 for numerical illustration, along with the RPV parameters determine neutrino 
masses and mixing angles( see section 2). The last two parameters are required to realize 
the top squark NLSP condition and the CCB condition respectively. It should be noted that 
the BR of the loop decay increases significantly for larger Ma- 
Top squark NLSPs having different masses are realized varying A t . For this set of RPC 
parameters the NLSP and the CCB conditions are satisfied for 940 < A t < 980. In addition 
we have used a common mass of L and R type sleptons Mj- = 350 and A T = 1000 for the 
computation of the BRs of the RPC modes. 

We then randomly generate 10 8 sets of the six RPV parameters in scenarios ( a\ and a 2 
) and filter out the ones allowed by the data in Table [1]. As discussed in section 3 in all 
allowed sets there is one dominant coupling ( A 233 or A 333 ). This hierarchy will obviously be 
reflected in the observed BRs. In Fig. 1 we present the BR of each of the three competing 
decay modes in Eq.(19) (including all leptonic channels), Eq.(22) and Eq.(23) ( including all 
possible / and /' combinations) and the corresponding number of allowed solutions out of 10 s 
randomly generated parameter sets. This figure clearly illustrates the competition among 
the three decay modes of the top squark NLSP. Here A t = 970GeV and = 181.5GeV. 
The RPV decay modes do not dominate although the combined BR of the three RPV modes 
is appreciable ( 20 % - 30 %) over most of the APS. It also follows from this figure that 
the number of solutions allowed by the v- data is indeed a tiny fraction of the total number 
of generated parameter sets. Thus the available oscillation data is already very resistive in 
spite of the large errors. 

In Ref.[9] top squark pair production followed by their RPV decays into e + e~ or fi + fi~ 
channels were studied for Tevatron RUN II. The minimum observable BR (MOBR) at Teva- 
tron for the decays into the e or fi channel have been estimated as a function of m^. For 
= 181.5GeV the MOBR of the (ti — > e + b) or (ti — > fi + b) channel is approximately 20 % 

It may be recalled that the analysis of [9] was conservative since only the leading order 
top squark pair production cross section and a total integrated luminosity of only 2000 pb^ 1 
were used. The next to leading order cross section is about 30% larger and accordingly a 
smaller estimate of the MOBR is expected. The total integrated luminosity during RUN II 
may be as large as 9000 pb~ l [23] which can further lower the MOBR. An improvement in 
the observability of RPV decays is also expected if ee, \i\x and e/x channels are simultaneously 
analyzed. 

For the solutions in which decays into the r lepton dominate, the combined BR of the 
modes involving lighter leptons are below the MOBR. Since this happens in a large region of 
the APS, it will be worthwhile to estimate the MOBR of this channel. On the other hand if 
the RPV decay into eor/i dominates, many allowed solutions have BRs close to the MOBR 
estimated in [9]. 
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For numerical illustrations we have considered the following points in the APS 

A) A' 133 = 3.72 x 10~ 5 , = 3.3 x 1(T 5 and A 333 = 8.5 x 10~ 5 ( scenario (en), A' 133 « A' 233 ) 

B) A' 133 = 4.2 x 10" 5 , A' 233 = 1.3 x 10 ~ 5 and A' 333 = 8.3 x 1(T 5 ( scenario (ai), A' 133 > A' 233 ) 

C) A' 133 = 4.2 x 10" 5 , A 233 = 8.6 x 1(T 5 and A' 333 = 1.8 x 10~ 5 ( scenario (a 2 ), A 233 > A' 333 ) 

The BRs of top squark decay modes into different leptonic channels for the parameter 
sets A) - C) are given in Table [4]. 



Decay modes 


ti -> e+b 


i"i -> fi+b 


ti -> T+b 


loop 


4 body 


BR(A) 


0.03 


0.024 


0.158 


0.69 


0.10 


BR(B) 


0.04 


0.004 


0.153 


0.70 


0.11 


BR(C) 


0.039 


0.163 


0.007 


0.69 


0.10 



Table 4: BRs of the competing decay modes of the lighter top squark in the gaugino domi- 
nated model ( Model 2, see text for the choice of parameters ) 

In Table [5] we present the number of events corresponding to different dileptonic final 
states at the Tevatron arising from RPV decays of both the produced top squarks. The 
number of events are computed for the BRs in Table [4], y^s = 2 TeV and an integrated 
luminosity of 9000 pb' 1 . For = 181.5 GeV the production cross section a(pp — > tit{) is 
0.41 pb as computed by CalcHEP v2. 1 [24]. The signal 

tJl -> l+bljb (24) 

is abbreviated as Lij, where i,j = 1,2,3 for e, p and r respectively. For i ^ j we considered 
dileptons of all charge combinations. 





Lu 


Ll2 


-^13 


L>22 


L23 


-^33 


# of events (A) 


3 


5 


34 


2 


27 


91 


# of events (B) 


5 


1 


44 





4 


86 


# of events (C) 


5 


46 


2 


97 


8 






Table 5: Typical sizes of opposite sign dileptons of different flavour combinations at the 
Tevatron from top squark pair production using the BRs in Table [4]. 

It is expected that the backgrounds can be suppressed to the desired level by employing 
standard kinematical cuts, b-tagging and by reconstructing the invariant masses of the two 
top squarks [9]. 

Another interesting signal arises if one of the produced top squarks decay via an RPV 
mode (a) while the other decays via the loop induced mode Eq. (22) followed by the LSP 
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decay leading to 

M* -> or (25) 

where i=l,2,3 as before and Lj denotes the above signal. The number of signal events for 
different % are presented in Table [6] using the same inputs as in Table [5]. We have included 
leptons of both signs in the signal. 





Li 


L 2 


L 3 


# of events (A) 


152 


121 


804 


# of events (B) 


206 


20 


789 


# of events (C) 


198 


829 


35 



Table 6: Typical sizes of various signals from top squark pair production at Tevatron RUN 
II when one of them decays into a RPV channel while the other into the loop induced mode 
followed by LSP decay. 

The top squark mass can be reconstructed by the invariant masses of the two hardest 
leptons and jets [9] in any Ly type event. The upper edge of the invariant mass spectrum 
of the two jets with lowest and next to lowest energy in any Lj type signal may provide 
information about the LSP mass. 

If the parameter K 2 , which sets the scale of the one loop mass matrix, is decreased keeping 
the parameters of the gaugino sector fixed, the allowed values of A^ 33 couplings increase and 
apparently larger BRs of the RPV modes are allowed. However, we shall illustrate the 
constrained nature of the model in the next paragraph and show that the above BRs cannot 
be arbitrarily large inspite of many free RPC parameters in the model. 

The parameter K 2 decreases for higher values of the common squark mass. In practice for 
a fixed set of gaugino parameters the common squark mass cannot be increased significantly 
without violating the top squark NLSP condition. Of course larger values of the A t parameter 
may restore the NLSP condition. But larger values of A t tend to violate the CCB condition. 
Finally the parameter M A can be increased to satisfy the CCB condition but as noted earlier 
that would enhance the loop decay width as well and the BRs of the RPV modes will still be 
suppressed. Thus once we know the parameters of the gaugino sector from complementary 
experiments and observe RPV decays of the top squark NLSP, the predicted BR of these 
modes cannot be made arbitrarily large by adjusting the common squark mass or Ma- 

The BR of the RPV decays increase significantly if we consider model 3 with mixed 
xf and x 2 - ^ n this case the parameter C increases substantially compared to it's typical 
magnitude in model 2 while the loop level mass matrix has a smaller K 2 due to a smaller 
higgsino mass parameter //. Thus the loop level mass matrix can be effective only for larger 
values of the trilinear RPV couplings. 

We demonstrate these effects with the same parameter set as above except that we take 
pL = 210.0 GeV. We present the corrsponding histogram in Fig. 2. It follows that the entire 
APS correspond to larger BRs of the RPV modes ( including all leptonic channels). 
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For illustrating the signals at the LHC we consider a new parameter space leading to a 
heavier top squark NLSP beyond the reach of the Tevatron. We have chosen m~ tl = 350.65 
GeV, Mi = 310, M 2 = 400, fi = 800, tan (3 = 5, M q = 550, Mj = 450, A b = 1000, A T = 1000, 
A t = 1350 and Ma = 500 where all masses and mass parameters are in GeV . At the LHC 
the top squark pair production cross section is 4.28 pb for this m^. 

We have chosen the following sets of RPV parameters consistent with the v data which 
reflects the same characteristics as in A) - C) listed above. 

A) A' 133 = 2.92 x 10~ 5 , A' 23 3 = 2.27 x 10~ 5 and A 333 = 6.27 x 10~ 5 

B) A' 133 = 3.4 x 10~ 5 , A 233 = 6.6 x 10~ 5 and A 333 = 1.83 x 10 ~ 5 

C) A' 133 = 3.39 x 10" 5 , A 233 = 1.92 x 10~ 5 andA 333 = 6.48 x 10~ 5 
The corresponding BRs are shown in Table [7]. 



Decay modes 


ti -> e + b 


h -> fi+b 


f"i -> T+b 


loop 


4 body 


BR(A) 


0.12 


0.075 


0.57 


0.23 


0.001 


BR(B) 


0.16 


0.59 


0.045 


0.21 


0.001 


BR(C) 


0.16 


0.05 


0.58 


0.21 


0.001 



Table 7: BRs of different decay modes of the lighter top squark in Model 2 ( see text for the 
choice of parameters ) 

Using the BRs in Table [7] and a representative integrated luminosity of 30/6 -1 one can 
easily calculate the number of various signal events. For example, we obtain 44696 L 22 and 
31817 L 2 events with the parameter set B. 

5. LSP decay 

As discussed in the introduction unless m^o > m t , LSP decay alone cannot provide detailed 
information about the underlying model of m v . 

Of course the decay of the LSP into the channel in Eq. (1) may provide circumstantial 
evidences in favour of an underlying RPV model of neutrino mass. For example, if Xi is 
assumed to be the LSP, then xtxl an d xtXi production followed by appropriate decay 
chains ending in LSP decays are indicative of an underlying model of neutrino mass [25]. In 
Ref. [25] the prospect of observing this signal at RUN II was studied. It was concluded that 
this signature can be probed up to mi/ 2 = 230 GeV(320 GeV) with an integrated luminosity 
of 2 fb" 1 ( 30 fb^ 1 ). Here m^j is the common gaugino mass at the GUT scale. However, 
the lepton number violating nature of the decay can not be established by the data. 

It may be noted that the pure RPC decay \X — > Xi^b, which may have a large BR if one 
of the bottom squark mass eigenstates happens to be lighter than the other squarks at large 
tan/5, has collider signatures very similar to the decay of Eq.(l). This is especially so if the 
LSP mass is much smaller than m^o which is quite common in models with non-universal 
gaugino masses. Thus one has to worry about the possibility of RPC SUSY faking the RPV 
signal. 
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Using model 1) (see section 4) we have chosen the following of RPC parameters : 
Mi = 710, M 2 = 800, fi = 395, tan (3 = 5, M~ q = 500, M~ x = 450, A b = 1000, A T = 1000, 
A t = 800 and M A = 300 where all masses and mass parameters are in GeV. This leads to 
m~ tx =391 GeV and m^o=388 GeV Since m^o > m t the following additional decay modes 
open up : 

(a) xi - ebt (26) 
(&)*?-> A* 6* (27) 
(c) x? - r 6 t (28) 

In Fig. 3 we present the BRs of the four competing decay modes of the LSP and the corre- 
sponding number of allowed solutions out of 10 8 sets of randomly generated RPV parameters. 

From the sets allowed by the neutrino data we have chosen the following RPV parameters 
: A' 133 = 5.43 x 10" 5 , = 10.19 x 10~ 5 A 333 = 0.87 x 10~ 5 . The resulting BRs are presented 
in column 2 of Table [8]. Here, the three neutrinos carry the fir- 



Decay modes 


BR (%)in Higgsino model 


BR (%)in Gaugino model 


Xi -> firbb 


12.5 


79.7 


Xi -»• teb 


19.3 


3.7 


Xi -»■ tfib 


67.7 


16.2 


Xi -> trb 


0.5 


0.4 



Table 8: BRs of different decay modes of the LSP in higgsino dominated model ( Model 1 ) 
and the gaugino dominated model (Model 2 ) ( see text for the choice of parameters). 

We next compare this result with the gaugino dominated scenario (Model 2). We choose 
Mi = 388, M 2 = 500, jj, = 800 and A t = 870 keeping all other RPC parameters same as 
that in the Higgsino type model. The choice of parameters is dictated by the fact that the 
masses of the LSP and different squarks remain practically the same in the two models being 
compared. The BRs in this model are presented in the last column of Table [8] The BRs 
of the modes involving charged leptons and the t quark are significantly enhanced in the 
Higgsino model because of the large top Yukawa coupling. 

We next study gaugino pair production followed by cascade decays with RPC and RPV 
parameters as in model(l) quoted above. We consider the following signal 

G a Gf3 -> ifhf^bttX (29) 

where G a Gp represents any pair of electroweak gauginos , i,j = 1,2,3 for e, \x and r respec- 
tively, Lij represents the number of this signal and X denotes any other particles produced. 
The number of signal events from different gaugino pair production at LHC are presented in 
Table [9]. The production cross-sections involved are calculated using CalcHEP v2.1 [24]. 
L12 includes e and \i events with all possible charge combinations and Ln/22 represents events 
with same sign di-leptons only. 

In Fig. 4 we present the LSP decay length vs the number of allowed solutions in Model 
1. It is seen that apart from a small region of the APS the decay will be inside the detector. 
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Gaugino pair 


a( pb) 


L n 


L>22 


L\2 


X? xl 


18.5 xl(T 3 


206 


2551 


5808 


xi xi 


68.2 xlO~ 3 


762 


9405 


21416 


X1X2 


10.2 xl(T 3 


113 


1406 


3200 


X1X2 


23.4 xlO~ 3 


261 


3226 


7348 


Xi Xi 


10.6 xlO" 3 


118 


1461 


3328 


v° V+ 
Al Al 


24.6 xlO~ 3 


274 


3392 


7724 



Table 9: Number of Ln,L 2 2 and L\ 2 events (see text) arising from various gaugino pair 
production at LHC followed by LSP decay 

6. Conclusions 

In this paper we have studied the APS of the RPV parameters in a model of neutrino mass 
subject to the constraints imposed by the neutrino oscillation data [14]. The model we have 
considered has three bilinear RPV couplings fii and three trilinear couplings A -33, where i 
is the lepton index, and we work in a basis where the sneutrino vev is zero. As expected 
from the upper bounds on the A -33 couplings obtained by the earlier analyses [12], we find 
for representative choice of the RPC parameters , the allowed magnitudes of these couplings 
are indeed very small (see section 3). Our analyses reveal that even the currently available 
//-oscillation data with relatively large errors are quite restrictive. Out of many randomly 
generated sets of RPV parameters consistent with the upperbounds on them [12], only a few 
are allowed by the data. 

Moreover, we also identify six generic scenarios, consistent with neutrino data, leading 
to distinctive collider signatures. These scenarios are listed as ai) - c 2 ) in section 3. In 
each scenario there is one small bilinear parameter /Zj and a characteristic hierarchy of A- 33 
couplings. Thus if a few direct RPV decay modes of any sparticle are observed their relative 
BRs would reflect the underlying model of neutrino mass. We have studied the decay modes 
of the top squark NLSP both at Tevatron RUN II and at the LHC ( section 4). Over the 
entire APS the BRs of the RPV decays are found to be significant for representative values 
of the RPC parameters and the hierarchy among them can potentially reveal the underlying 
model of neutrino mass. The LSP decays can provide similar information if the LSP mass 
is larger than m t (section 5). 

Hopefully RPV decays in the above channels will be seen at hadron colliders and at least 
some parameters of the RPC sector can be measured kinematically ( some examples are 
discussed in section 4 ). The accurate measurement of the relative BRs at the ILC will then 
provide an exciting program for probing the origin of neutrino mass. 

In this paper we have chosen a particular set of RPV parameters in the physical basis at 
the weak scale for explaining the current v- oscillation data. Following the standard practice 
we have further assumed that the magnitudes of all other RPV parameters are negligible. 
There is, however, another exciting possibility. Let a few relatively large RPV parameters 
not directly related to the v - sector, be generated at a high scale (say M GUT ) by any suitable 
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mechanism ( for examples by effective operators [26] ). The renormalization group ( RG ) 
evolution of these parameters down to the weak scale can then generate the parameters ^ 
and A^ 33 having magnitudes severely suppressed compared to the input parameters [?]. This 
happens due to the flavour violation ( non-diagonal Yukawa couplings matrices ) inevitably 
present in the quark sector. It was shown in [27, 28] that any three input couplings at Mqvt 
with different lepton indices chosen from A- 13 and A- 23 , i=l,2,3, can induce a viable v- mass 
matrix at the weak scale. Thus at the weak scale there may be some relatively large A' 
couplings in addition to the small parameters underlying v- physics. Moreover when the 
quark mass matrices are diagonalized at the weak scale rotations on the quark fields may 
further generate new RPV couplings [28] relevant for v- physics from the above relatively 
large A'- type couplings. Of Course such rotations may also induce other A' -type couplings 
leading to unacceptable flavour changing neutral currents [29]. Care should, therefore, be 
taken in choosing the input parameters at M GUT . Because of the relatively large A'- type 
couplings the RPV collider phenomenology and some rare decays of K-mesons and r-leptons 
would be rather spectacular in these models. Several specific models and their associated 
phenomenologies were discussed with numerical illustrations in [28]. 
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Figure 1: The branching ratios(%) of the three competing decay modes vs the number of allowed 
solutions in Model 2 ( see text for the parameters used ). 
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Figure 2: The branching ratios(%) of the two competing decay modes vs the number of allowed 
solutions in Model 3 ( see text for the parameters used ). 
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Figure 3: The branching ratios (%) of competing LSP decay modes vs the number of allowed 
solutions in Model 1 (see the text for the parameters used ). 
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Figure 4: Neutralino decay length vs the number of allowed solutions in Model 1 (see the text for 
the parameters used ). 
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